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I. INTRODUCTION 

Synchronisation of nonlinear oscillators to irregular external signals is an interesting problem of importance in 
physics, biology, applied science, and engineering [ll4l0|. The key difference to synchronisation to a periodic external 
signal is the lack of a simple functional relationship between the input signal and the synchronised output signal, 
making the phenomenon much less evident [H. [28j. Rather, synchronisation is detected when two or more identical 
uncoupled oscillators driven by the same external signal but starting at different initial states have identical long-term 
responses. This is equivalent to obtaining reproducible long-term response from a single oscillator driven repeatedly 
; i ■ by the same external signal, each time starting at a different initial state. Hence, synchronisation to irregular external 
signals is also known as reliability Q or consistency 0, and represents the ability to encode irregular signals in a 
reproducible manner. 

Recent studies have shown that nonlinear oscillators can exhibit interesting responses to stochastic external signals. 
Typically, a small amount of external noise causes synchronisation [l|, Ch.7], jl, Hj]. However, as the strength of 
external noise increases, there can be a loss of synchrony in oscillators with amplitude-phase coupling (also known 
as shear, nonisochronicity or amplitude-dependent frequency) 0, Ch.7], [il-fiol fl3 - fl^ | . Mathematically, loss of noise 
synchrony, consistency or reliability is a manifestation of a stochastic bifurcation of a random attractor. 

This chapter gives a definition of noise synchronisation in terms of random pullback attractors and studies 
synchronisation-desynchronisation transitions as purely noise-induced stochastic bifurcations. This is in contrast 
to the effects described in Chapter 2, where noise is used to control or regulate the dynamics that is already present 
in the noise-free system. We focus on a single-mode class-B laser model and the Landau-Stuart model (Hopf normal 
form with shear In section [VTl numerical analysis of the locus of the stochastic bifurcation in a three-dimensional 
l— "™ ', parameter space of the 'distance' from Hopf bifurcation, amount of amplitude-phase coupling, and external signal 
strength reveals a simple power law for the Landau-Stuart model but quite different behaviour for the laser model. 
£>. \ In section IVII| the analysis of the shear-induced stretch-and-fold action that creates horseshoes gives an intuitive 
£SJ ■ explanation for the observed loss of synchrony and for the deviation from the simple power law in the laser model. 
Experimentally, stochastic external forcing can be realised by optically injecting noisy light into a (semiconductor) 
laser as described in section HVl While bifurcations of random pullback attractors and the associated synchronisation- 
desynchronisation transitions have been studied theoretically, single-mode semiconductor lasers emerge as interesting 
, candidates for experimental testing of these phenomena. 
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II. CLASS-B LASER MODEL AND LANDAU-STUART MODEL 

A class-B single- mode laser [l5[ without noise can be modelled by the rate equations [lij : 

dE 

— = iAE + g-f(l - ia)NE, (1) 
at 

dN . , „, 

— = J-N-(l + gN)\E\ 2 , (2) 

which define a three-dimensional dynamical system with a normalised electric field amplitude, E € C, and normalised 
deviation from the threshold population inversion, JVgl, such that N — — 1 corresponds to zero population inversion. 
Parameter J is the normalised deviation from the threshold pump rate such that J = — 1 corresponds to zero pump 
rate. The linewidth enhancement factor, a, quantifies the amount of amplitude-phase coupling, A is the normalised 
detuning (difference) between some conveniently chosen reference frequency and the natural laser frequency, 7 = 500 
is the normalised decay rate, and g = 2.765 is the normalised gain coefficient (l6| . 

System (JTH2]) is § 1 -equivariant, meaning that it has rotational symmetry corresponding to a phase shift E —> Ee 1 ^, 
where < <f> < 2ir. For JgR, there is an equilibrium at (E, N) — (0, J) which represents the "off" state of the laser. 
This equilibrium is globally stable if J < and unstable if J > 0. At J = 0, there is a Hopf (A 7^ 0) or pitchfork 
(A = 0) bifurcation defining the laser threshold. Moreover, if J > 0, the system has a stable group orbit in the form 
of periodic orbit for A ^ or a circle of infinitely many non-hyperbolic (neutrally stable) equilibria for A = 0. In 
this paper, we refer to this circular attractor as the limit cycle. The limit cycle is given by (|-E| 2 ,iV) = ( J, 0) and 
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represents the "on" state of the laser. Owing to the S 1 -symmetry, the Floquet exponents of the limit cycle can be 
calculated analytically as eigenvalues of one of the non-hyperbolic equilibria for A = 0. Specifically, if 



< J < 



4 7 (l - y/1 - l/(2 7 )) 



9 x 10" 



the overdamped limit cycle has three real Floquet exponents 

Mi = 0, a*2,3 =-a±b, 

and if 



(3) 



4 7 (l - VI - V(27)) 



1 4 7 (l + y/1 - l/(2 7 )) 
— < J < * '- 
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9 9 
the underdamped limit cycle has one real and two complex-conjugate Floquet exponents 

Mi = 0, M2,3 = —a± ib, 

where 



1446, 



(4) 



and b = y/\a 2 - 2gj J\ > 0. 



In the laser literature, the decaying oscillations found for pump rate in the realistic range J £ (9 x 10 _5 ,20) are 
called relaxation oscillations. (This should not be confused with a different phenomenon of self-sustained, slow-fast 
oscillations.) Finally, even though the laser model (Q][2]) is three dimensional, it cannot admit chaotic solutions due 
to restrictions imposed by the rotational symmetry. 

Using centre manifold theory [17| . the dynamics of (jT] [SJ) near the Hopf bifurcation can be approximated by the 
two-dimensional invariant centre manifold 



W c = {(E, N) £ M 3 : N = J - \E\ 2 }, 



on which ((TH21) reduces to 

1 dE 

91 dt 

After rescaling time and detuning, 

we obtain the Landau-Stuart model 

dE 
dt 



J + i ( A _ a (J - \E\ 2 ) 



E- E\E[ 



t = tgi and A = A/ (# 7 ), 



J + i A-a(J-|£| ) 



E-E\E\ 



(5) 



that is identical to the Hopf normal form }17| except for the higher order term, ia(J — \E\ 2 )E, representing amplitude- 
phase coupling. Since this term does not affect stability properties of (0, it does not appear in the Hopf normal form. 
However, in the presence of an external forcing, f ex t(t), this term has to be included because it gives rise to qualitatively 
different dynamics for different values of a. If J > 0, the Landau-Stuart model has a stable limit cycle with two 
Floquet exponents 



fj,X = Q and fj,2 = — 2J. 



(6) 



III. THE LINEWIDTH ENHANCEMENT FACTOR AND SHEAR 



The linewidth enhancement factor, a, quantifying the amount of amplitude-phase coupling for the complex-valued 
electric field, E, is absolutely crucial to our analysis. Its physical origin is the dependence of the semiconductor 
refractive index, and hence the laser-cavity resonant frequency, on the population inversion [l5l [l8j . A change in the 
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(a) a = 




-2 -2 Mm 



FIG. 1: (Black) The limit cycle representing the "on" state of the laser for J = 1 and (gray) isochrones for three different 
points on the limit cycle as defined by Eq. (0 for (a) a = and (b) a — 2. 
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FIG. 2: An experimental setup for detecting noise synchronisation in lasers. 
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electric field intensity, <5|-E| 2 , induces a change, SN, in population inversion [Eq. ((2J)]. The resulting change in the 
refractive index shifts the cavity resonant frequency. The ultimate result is a change of —ag^/SN in the instantaneous 
frequency of the electric field denned as d(&Tg(E))/dt. 

Mathematically, amplitude-phase coupling is best illustrated by an invariant set associated with each point, q, on 
the limit cycle. For a point g(0) on a stable limit cycle in a n-dimensional system, this set is denned as 

{x(0) G R" : x(t) -» q(t) as t -> oo}, (7) 

In the laser model (JTH2]) and Landau-Stuart model ((2]), isochrones are logarithmic 



and is called an isochrone 
spirals that satisfy 



To see this, define a phase 



arg(£) + a In \E\ = C, where C G (0, 2tt]. 



# = arg(_E) + a In \E\ 



(8) 



(9) 



and check that d^/dt is constant and equal to A for (jT|-[2]) and A for ||5j). This means that trajectories for different 
initial conditions with identical initial phase, ^(0), will retain identical phase, ^(i), for all time t. Since the limit 
cycle is stable, all such trajectories will converge to the limit cycle, where they have the same |i?(t)|. Then, Eq. © 
implicates that all such trajectories have the same axg(E(t)) and hence converge to just one special trajectory along 
the limit cycle as required by (J7J. 

Isochrones of three different points on the laser limit cycle are shown in Fig.[TJ Isochrone inclination to the direction 
normal to the limit cycle at q(0) indicates the strength of phase space stretching along the limit cycle. If a = 0, 
trajectories with different \E\ > rotate around the origin of the i?-plane with the same angular frequency giving 
no isochrone inclination and hence no phase space stretching [Fig.QJa)]. However, if \a\ > 0, trajectories with larger 
\E\ rotate with higher angular frequency giving rise to isochrone inclination and phase space stretching [Fig. H£b)]. 
Henceforth, we refer to amplitude-phase coupling as shear. 



IV. DETECTION OF NOISE SYNCHRONISATION 



There are at least two approaches to detecting synchronisation of a semiconductor laser to an irregular external 
signal. One approach involves a comparison of the responses of two or more identical and uncoupled lasers that 
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are driven by the same external signal. The other approach involves a comparison of the responses of a single laser 
driven repeatedly by the same external signal Here, we consider responses of M uncoupled lasers with intrinsic 
spontaneous emission noise that are subjected to common optical external forcing, f e xt(t), fl6L [20j| : 

dE ■ 

= iAE j +gy(l-ia)N J E j + fEj(t) + f aBt (t) (10) 

^ = J-Ni-il+gN^Etf + f^it), (11) 
3 = 1,2,. ..,M, 

The lasers are identical except for the intrinsic spontaneous emission noise that is represented by random Gaussian 
processes 

fE j (t) = f§ j (t) + if I Ej (t) and f NJ (t), 
that have zero mean and are delta correlated 

(/**(*)) = (fNj(t)) = 0, 

(/£(*)/£#)) = (/ii(*)/ii(0> = ^6(1 - 0, (12) 
(fm(t)f Nj (t')) = 2D N S. lJ S(t-t'). 

Here, #jj is the Kronecker delta and <5(t — t') is the Dirac delta function. In the calculations we use De = 0.05 and 
D N = 3.5 x 10" 8 QJ]. 

To measure the quality of synchronisation we introduce the order parameter, Im (t) , and the average order parameter, 
(I M ), as 



(7 M ) = lim )= ( I M (t) dt = lim ± / 



M 



2 

(13) 



The physical meaning of Im (t) and (Im) is illustrated in Fig.[5J If M identical lasers are placed at an equal distance 
from a small (the order of a wavelength) spot and their light is focused onto this spot, then Jjut (t) and (Im) are 
the instantaneous and average light intensity at the spot, respectively. A single laser oscillates with a random phase 
owing to spontaneous emission noise so that, for independent lasers, (Im) is proportional to M times the average 
intensity of a single laser. This follows directly from Eq. (1131) assuming lasers with identical amplitudes, |.Ej(t)|, and 
uncorrelated random phases, axg(Ej(t)), However, when the lasers oscillate in phase, one expects (Im) to be equal 
M 2 times the average intensity of a single laser. This follows directly from Eq. (fT3|) assuming lasers with identical 
amplitudes and phases. We speak of synchronisation when (Im) ~ M 2 , different degrees of partial synchronisation 
when M < (Im) < M 2 , and lack of synchronisation when (Im) ~ M. Note that (Im) > M 2 indicates trivial 
synchronisation, where the external forcing term, f ex t{t) becomes 'larger' than the oscillator terms on the right-hand 
side of Eq. (fl~0|) . For comparability reasons, we now briefly review the case of a monochromatic forcing and then move 
on to the case of stochastic forcing. 

Let us consider a monochromatic external forcing 

f ext (t) = Ke™^\ 

where K G M is the forcing strength and v ext is the detuning (difference) between the reference frequency chosen 
for A in Eq. ([1]) and the forcing frequency. Such an external forcing breaks the S^-symmetry and can force each 
laser to fluctuate in the vicinity of the well-defined external forcing phase, v ext t, as opposed to a random walk. This 
phenomenon was studied in [2 11 ] as a thermodynamic phase transition. Figure [3fa) shows (Im) versus K, for an 
external forcing resonant with the laser, 



Because of the intrinsic spontaneous emission noise, the forcing amplitude has to reach a certain threshold before 
synchronisation occurs. For a — 0, a sharp onset of synchronisation at K m 10 -3 is followed by a wide range of K 
with synchronous behaviour where (Im) = M 2 (If r ). Here, (If r ) is the average intensity of a single laser without 
forcing. At around K = 10 2 , (Im) starts increasing above M 2 (If r ). Whereas lasers still remain synchronised, this 
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FIG. 3: The average order parameter as defined by Eq. (|13|) for M — 50 uncoupled lasers with common (a) monochromatic 
and (b) white noise external forcing vs. the forcing strength for (dashed) a = and (solid) a — 3; J = 5 and v ext — A in panel 
(a). Adapted from [13] with permission. 
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FIG. 4: The shape of the probability distribution of 7m (t) for white noise external forcing; (dashed) a = and (solid) a = 3. 
From (a) to (c) D ext = 10°, 10 1 , and 10 3 . Adapted from pj with permission. 

increase indicates that the external forcing is no longer 'weak'. Rather, it becomes strong enough to cause an increase 
in the average intensity of each individual laser. A very different scenario is observed for a — 3. There, the onset of 
synchronisation is followed by an almost complete loss of synchrony just before (7m) increases above M 2 (If r ). The 
loss of synchrony is caused by externally induced bifurcations and ensuin g ch aotic dynamics. These bifurcations have 
been studied in detail, both theoretically [l5|, l22l - [24j and experimentally |25| | , and are well understood. 

The focus of this work is synchronisation to white noise external forcing represented by the complex random process 
that is Gaussian, has zero mean, and is delta correlated 

/«*(*) = /£*(*)+</«*(*). 

(/«t(t)) = (f* t (t)fLM = o, (14) 

</£*(*)/£*(*')) = (fLt(t)fLt(t'))=D ext 5(t-t>). 

White noise synchronisation is demonstrated in Fig. [3jb) where we plot (7m) versus D ext . For a = 0, a clear onset 
of synchronisation at around D ext = 10~ 3 is followed by synchronous behaviour at larger D ext . In particular, there 
exists a range of D ext where white noise external forcing is strong enough to synchronise phases of intrinsically noisy 
lasers but weak enough so that each individual laser has small intensity fluctuations and its average intensity remains 
unchanged. In the probability distributions for Jm(£) hi Fig- IHa _ b), the distinct peak at 7m(0 ~ M 2 (Ifr) and 
a noticeable tail at smaller I]\i(t) indicate synchronisation that is not perfect. Rather, synchronous behaviour is 
occasionally interrupted with short intervals of asynchronous behaviour owing to different intrinsic noise within each 
laser. For D ext > 10 2 the external forcing is no longer 'weak' and causes an increase in the intensity fluctuations and 
the average intensity of each individual laser. Although the lasers remain in synchrony, (7m) increases above M 2 (If r ) 
[Fig-EIb)] and exhibits large fluctuations [Fig. H[c)] as in the asynchronous case. A very different scenario is observed 
again for a = 3. There, the onset of synchronisation is followed by a significant loss of synchrony for D ext £ (4, 100). 
In this range of the forcing strength, one finds qualitatively different dynamics for a = and a — 3 as revealed by 
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different probability distributions in Fig. [4] (b) . 

Interestingly, comparison between (a) and (b) in Fig. [3] shows that some general aspects of synchronisation to a 
monochromatic and white noise external forcing are strikingly similar. In both cases there is a clear onset of synchro- 
nisation followed by a significant loss of synchrony for sufficiently large a, and subsequent revival of synchronisation 
for stronger external forcing. However, the dynamical mechanism responsible for the loss of synchrony in the case of 
white noise external forcing has not been fully explored. 



V. DEFINITION OF NOISE SYNCHRONISATION 



The previous section motivates further research to reveal the dynamical mechanism responsible for the loss of syn- 
chrony observed in Fig. EJb). To facilitate the analysis, we define synchronisation to irregular external forcing within 
the framework of random dynamical systems. Let us consider a n-dimensional, nonlinear, dissipative, autonomous 
dynamical system, referred to as unforced system 

^ =/(x,p), (15) 

where x £ W l is the state vector and p € M. k is the parameter vector that does not change in time. An external forcing 
is denoted with f ex t(t), and the corresponding non-autonomous forced system reads 

Y t =/(z,P)+/e*t(*)- (16) 

Let x(t, to, Xq) denote a trajectory or solution of (|16[) that passes through xq at some initial time to- In situations 
where explicitly displaying the initial condition is not important we denote the trajectory simply as x(t). For an 
infinitesimal displacement Sx(0) from x(0, to, xq), the largest Lyapunov exponent along x(t,to,xo) is given by 

A™, = km i In |^T- (17) 
t^oo t \Sx(0)\ 

If the external forcing, fext(t)Jis stochastic, Eq. ()16[) defines a random dynamical system where X max does not depend 
on the noise realisation, f ext [26| . Furthermore, we define: 

Definition 1. An (self-sustained) oscillator is an unforced system (|15|) with a stable hyperbolic limit cycle. 

Definition 2. An attractor for the forced system (fTB| with stochastic forcing f e xt(t) is called a random sink (rs) if 
X m ax < 0, and a random strange attractor (rsa) if X max > 0. 

Definition 3. A stochastic d-bifurcation is a qualitative change in the random attractor when X ma x crosses through 
zero [H, Ch.9]. 

Definition 4. An oscillator is synchronised to stochastic forcing f e xt(t) on a bounded subset D c R™ if the corre- 
sponding forced system (I16|) has a random sink in the form of a unique attracting trajectory, a(t,f ex t), such that 



lim \x(t,t ,x ) - a(t, f ext )\ ^> 0, 
to — >— oo 



for fixed t > to and all xo G D. 



By Definition 1, an unforced oscillator has zero X ma x on an open set of parameters. In the presence of stochastic 
external forcing, X max becomes either positive or negative for typical parameter values [13, [26| and remains zero only 
at some special parameter values defining stochastic d-hifur cations. Synchronisation in Definition 4 is closely related 
to generalised synchronisation [28l.[2^| or weak synchronisation [30l | — a phenomenon that requires a time-independent 
functional relationship between the measured properties of the forcing and the oscillator [3l|. Following Refs. [2(1 
Ch.9] and [32j . we used in Definition 4 the notion of pullback convergence where the asymptotic behaviour is studied 
for to — > — oo and fixed t. (For a study of different notions of convergence in random dynamical systems we refer 
the reader to Ref. [27].) Whereas X max does not depend on the noise realisation, f ex t, random sinks and random 
strange attractors do depend on f ext . Hence the f ex t dependence in a{t, f ex t) in Definition 4. Since X max < does not 
imply a unique attracting trajectory it is not sufficient to show synchronisation as defined in Definition 4. In gen eral, 
there can be a number of coexisting (locally) attracting trajectories that belong to a global pullback attractor [321 ] . In 
such cases, one can choose D to lie in the basin of attraction of one of the locally attracting trajectory and speak of 
synchronisation on D. 
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FIG. 5: Bifurcation diagram for the white noise forced laser model (|18H19[) in the (D ex t, J)-plane for (a) a = and (b) a = 3. 
rs stands for random sink, rsa stands for random strange attractor, the two black curves d denote stochastic d-bifurcation, and 
colour coding is for \ ma x- An unforced laser has a Hopf bifurcation at {D ex t = 0, J = 0). Adapted from [l2T ] with permission. 



VI. SYNCHRONISATION TRANSITIONS VIA STOCHASTIC D-BIFURCATION 



To facilitate the analysis we make use of Definitions 1-3 and, henceforth, consider noise synchronisation in the laser 
model with white noise external forcing 



dE 

~dt 

dN 

~dt 



iAE + 57 (1 - ia)NE + f ext (t), 
J-N-(l + gN)\E\ 2 , 



(18) 
(19) 



but without the intrinsic spontaneous emission noise. Now, owing to the absence of intrinsic noise, Definition 4 
is equivalent to the synchronisation detection scheme chosen in Section IIVI More specifically, the evolution of M 
trajectories starting at different initial conditions for a single laser with external forcing is the same as the evolution 
of an ensemble of M identical uncoupled lasers with the same forcing, where each laser starts at a different initial 
condition. A random sink for (|18H19J) in the form of a unique attracting trajectory, a(t,f ex t), makes trajectories 
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Re(E) Re(E) Re(E) 

FIG. 6: Pullback convergence to (a-c) a random sink for D^ x t = 0.1 and (d-f) a random strange attractor for D ext = 0.5 in 
the white noise forced laser model (|18H19|I with a = 3 and J = 1. Shown are snapshots of 10 000 trajectories in projection onto 
the complex _E-plane at time t = 30. The initial conditions are uniformly distributed on the _B-plane with N = at different 
initial times to=(a,d) 29, (b,e) 28, and (c,f) 0. 



for different initial conditions converge to each other. In an ensemble of M identical uncoupled lasers with common 
forcing this means that (Im) ~^ M 2 (If r ) in time so that synchronisation is detected. A random strange attractor 
for (IT8T[TI))) , where nearby trajectories separate exponentially fast because \ max > 0, implies (Im) < M 2 (If r ) so that 
incomplete synchronisation or lack of synchronisation is detected. 

Figure [5] shows effects of white noise external forcing on the sign of the otherwise zero A m(M hi an unforced laser. 
For a = 0, external forcing always shifts Xmax to negative values meaning that the system has a random sink for 
D ex t > and J > [Fig.JSfa)]. Additionally, this random sink is a unique trajectory, a(t, fext), meaning that the laser 
is synchronised to white noise external forcing. However, for a = 3 there are two curves of stochastic d-bifurcation 
where X max crosses through zero [Fig. [SJb)]. Noise synchronisation is lost for parameter settings between these two 
curves, where X ma x > indicates a random strange attractor . Pullback convergence to two qualitatively different 
random attractors found for a = 3 is shown in Fig. [6l At fixed time t — 30, we take snapshots of 10 000 trajectories 
for a grid of initial conditions with different initial times t$. In (a-c), trajectories converge in the pullback sense to 
a random sink. The random sink appears in the snapshots as a single dot whose position is different for different t 
or different noise realisations fext- I n (d-f), trajectories converge in the pullback sense to a random strange attractor 
that appears in the snapshots as a fractal-like structure. This structure remains fractal-like but is different for different 
t or different noise realisations fext- 



A. Class-B Laser Model vs. Landau-Stuart Equations 

The stochastic d-bifurcation uncovered in the previous section has been reported in biological systems [81-fToL [l3l llil ] 
and should appear in a general class of oscillators with stochastic forcing. Here, we use the laser model in conjunction 
with the Landau-Stuart model to address its dependence on the three parameters: D ex t, J, and a, and to uncover its 
universal properties. With an exception of certain approximations [10(, this problem is beyond the reach of analytical 
techniques and so numerical analysis is the tool of choice. 

To help identifying effects characteristic to the more complicated laser model, we first consider the Landau-Stuart 
model with white noise external forcing 

E - E\E\ 2 + f ext (i), (20) 
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FIG. 7: A three-parameter study of stochastic d-bifurcation (a-b) in the white noise forced Landau-Stuart model ((20} and 
(c-d) the white noise forced laser model (|18H19[) . To facilitate the comparison between (I20|l and (|18II19[) . we plotted results 
for (fTMl9)> with the rescaled forcing strength, y /, 2D ext /g'y. Adapted from [l^] with permission. 



where, for the rescaled time t, the external forcing correlations become 

(fUt)fe Xt (t')) = (fUt)fLt(t')) = —s(t-t>). 

91 

Figure [7j a-b) shows the dependence of the d-bifurcation on J, D ext , and a in Eq. (|20|) . In the three-dimensional 
( J, D ext , a)-parameter space, the two-dimensional surface of d-bifurcation appears to originate from the half line 
{D ex t — 0, J = 0,O! > 5.3) of the deterministic Hopf bifurcation, has a ridge at a m i n ~ 5.3, and is asymptotic to 
a m 9 with increasing D ext [Fig. [7(b)]. Furthermore, numerical results in Fig. [JJb) suggest that the shape of the 
d-bifurcation curve in the two-dimensional section (D extl a) is independent of J. As a consequence, for fixed a within 
the range a £ (5.3,9) one finds two d-bifurcation curves in the (D ext , J)-plane [Fig. [7(a)] that are parametrised by 

Jj = C J (a) v / 2D ext , where j = 1, 2, (21) 

and bound the region with a random strange attractor. Since C\{a m i n ) = Ci{a m in) — 1, these two curves merge 
into a single curve 

J = ^2D ext , (22) 

when a — ct m i n . On the one hand, for a < a m i n , the region with a random strange attractor disappears from the 
{D ext , J)-plane. On the other hand, for a > 9, there is just one d-bifurcation curve in the (D ext , J)-plane, meaning 
that the region with a random strange attractor becomes unbounded towards increasing D ext [Fig. [7(b)]. 

Similar results are expected for any white noise forced oscillator with shear that is near a Hopf bifurcation, and for 
'weak' forcing. This claim is supported with numerical analysis of the laser model flT51HT)|) in Fig. [7( c-d). For a fixed 
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FIG. 8: J versus (solid) the two non-zero Lyapunov exponents, Re(/i2) and Re(/i3), for the limit cycle in the laser model HHU, 
and (dashed) the nonzero Lyapunov exponent, /X2, for the limit cycle in the Landau-Stuart model ((5|. Adapted from with 
permission. 

a, Eq. (|20l) and Eqs. ([TBHT9")) give identical results if the forcing is weak enough but significant discrepancies arise with 
increasing forcing strength. First of all, it is possible to have one-dimensional sections of the (D ext , J)-plane for fixed 
J with multiple uplifts of X max to positive values [black dots for J < 10~ 6 in Fig.[7£c)]. Secondly, the parameter region 
with a random strange attractor for Eqs. (fl"8Hl"9|) expands towards much lower values of a > a m i n ~ 1 [Fig. E(d)]. 
Thirdly, the shape of the two-dimensional surface of d-bifurcation in the laser model becomes dependent on J and 
has a minimum rather than a ridge. As a consequence, although the stochastic d-bifurcation seems to originate from 
the half line (D ext = 0, J = 0, a > 5.3), it will appear in the (D ext , J)-plane even for a £ (1,5.3) as a closed and 
isolated curve away from the origin of this plane [Fig. EJc)]. Finally, the region of random strange attractor remains 
bounded in the (D ext , J)-plane even for large a. 

To unveil the link between the transient dynamics of unforced systems and the forcing-induced stochastic d- 
bifurcation, we plot J versus Lyapunov exponents in Fig. [SJ note that Lyapunov exponents, A^, and Floquet ex- 
ponents, Hi, are related by A,; = Re[/ij]. A comparison between Figs. [7] and [8] shows strong correlation between the 
relaxation towards the limit cycle and the d-bifurcation. In the Landau-Stuart model (JSJ, the linear relation © 
between J and the non-zero Lyapunov exponent, fi2 (dashed line in Fig. [8|, results in a linear parametrisation (|2ip 
of d-bifurcation curves in the (J, \ / 2D ext )-p\a,ne [Fig. 0[a)]. In the laser model (HH2]), the nonlinear relation 
between J and the non-zero Lyapunov exponents, Re(/i2) and Ke(fi^) (solid curves in Fig. |5}, results in a very similar 
nonlinear parametrisation of d-bifurcation curves in the (J, v / 2£) e2;t )-plane [Fig. 0c)]. The splitting up of the chaotic 
region bounded by the black dots for J < 9 x 10 5 in Fig. [2c) is related to two different eigendirections normal to the 
limit cycle with significantly different timescales of transient dynamics towards the limit cycle (the two corresponding 
Lyapunov exponents are shown in solid in Fig. [5J . Finally, the appearance of relaxation oscillations in the laser system 
is associated with a noticeable expansion of the chaotic region, in particular, towards smaller a. 

VII. NOISE-INDUCED STRANGE ATTRACTORS 

Complicated invariant sets, such as strange attractors, require a balanced interplay between phase space expansion 
and contraction [l7| . If phase space expansion in certain directions is properly compensated by phase space contraction 
in some other directions, nearby trajectories can separate exponentially fast (X max > 0) and yet remain within a 
bounded subset of the phase space. 

It has been recently proven that, when^ suitably perturbed, any stable hyperbolic limit cycle can be turned into 
'observable' chaos (a strange attractor) [33]. This result is derived for periodic discrete-time perturbations (kicks) 
that deform the stable limit cycle of the unkicked system. The key concept is the creation of horseshoes via a stretch- 
and-fold action due to an interplay between the kicks and the local geometry of the phase space. Depending on the 
degree of shear, quite different kicks are required to create a stretch-and-fold action and horseshoes. Intuitively, it 
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FIG. 9: Time evolution of sets of initial conditions showing the creation of horseshoes in the phase space of a suitably kicked 
laser model with no shear (a = 0) . The sets of initial conditions arc (a) the stable circle and (d) boxes containing parts of the 
circle. Shown are phase portraits [(a) and (d)] before, [(b) and (e)] immediately after the first kick, and [(c) and (f)] some time 
after the first kick. 

can be described as follows. In systems without shear, where points in phase space rotate with the same angular 
frequency about the origin of the complex .E-plane independent of their distance from the origin, the kick alone has 
to create the stretch-and-fold action. This is demonstrated in Fig. [9j Horseshoes are formed as the system is suitably 
kicked in both radial and angular directions and then relaxes back to the attractor [the circle in Fig. Eta)] of the 
unkicked system. Repeating this process reveals chaotic invariant sets. However, showing rigorously that a specific 
kick results in 'observable' chaos is a non-trivial task [33[. In the presence of shear, where points in phase space 
rotate with different angular frequencies depending on their distance from the origin, the kick does not have to be so 
specific or carefully chosen. In fact, it may be sufficient to kick non-uniformly in the radial direction alone, and rely 
on natural forces of shear to provide the necessary stretch-and-fold action. 

These effects are illustrated in Fig. [TU] for the single laser model (p]-[3]) with non- uniform kicks in the radial direction 
alone for a = (no shear) and a = 2 (shear). There, we set A = and refer to the stable limit cycle [dashed 
circle in Fig. ITUT a)] as T. Kicks modify the electric field amplitude, \E\, by a factor of 0.8sin[4 arg(_E)] at times 
t = 0, 0.25, 0.5, and 0.75, but leave the phase, axg(E), unchanged. For a = each point on the black curve spirals 
onto r in time but remains within the same isochrone defined by a constant electric-field phase, arg(_E) = arg(_E(0)). 
Hence, the black curve does not have any folds at any time. However, for a = 2, a kick moves most points on the 
black curve to different isochrones so that points with larger amplitudes \E\ rotate with larger angular frequencies. 
This gives rise to an intricate stretch-and-fold action. Folds and horseshoes can be formed under the evolution of the 
flow even though the kicks are in the radial direction alone. 

In the laser model, stretch-and-fold action is significantly enhanced by the spiralling transient motion about T. For 
J > 10 _1 , the laser model (JTH^J) and the Landau-Stuart model (JS|) have nearly identical relaxation timescales toward 
r (Fig. [5]). However, owing to one additional degree of freedom and oscillatory relaxation the instantaneous 
stretching along T in the three-dimensional laser vector field (Q~H2]) can be much stronger compared to the planar 
vector field ([5]), especially at short times after the perturbation. This effect is illustrated in Fig. [TT] by the time 
evolution of the phase difference, arg(_Ei(i)) — arg(£ , 2(i)), between two trajectories, 1 and 2, starting at different 
isochrones for a — 3. Since both vector fields have identically shaped isochrones ([5]), the phase difference converges 
to the same value as time tends to infinity. However, at small t, the oscillatory phase difference for (Q~H2]) can exceed 
significantly the monotonically increasing phase difference for ([5]) (compare solid and dotted curves in Fig. Ill[) . 

It is important to note that the rigorous results for turning stable limit cycles into chaotic attractors are derived for 
periodic discrete-time perturbations. Stochastic forcing is a continuous-time perturbation, meaning that the analysis 
in [33| | cannot be directly applied to our problem. Nonetheless, such analysis gives a valuable insight as to why random 
chaotic attractors appear for a sufficiently large, and it helps distinguish effects of stochastic forcing. 

Here, we demonstrated that purely additive white noise forcing is sufficient to induce random strange attractors 
in limit cycle oscillators. Furthermore, numerical analysis in Sec. I VI Al shows that, in the case of stochastic forcing, 
creation of strange attractors requires a different balance between the amount of shear, relaxation rate toward the 
limit cycle, and forcing strength, as compared to periodic forcing. Unlike in the case of discrete-time periodic forcing, 
the shear has to be strong enough, |a| > C > 0, to allow sufficient stretch-and-fold action. Provided that the shear is 
strong enough, the stochastic forcing strength needs to be at least comparable to the relaxation rate toward the limit 
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a = (a) a = 2 a = (b) a = 2 




a = (c) a = 2 a = G (d) a = 2 




FIG. 10: Snapshots at times (a) t = 0, (b) i = 0.35, (c) t = 0.8, and (d) t = 1 showing the time evolution of 15000 trajectories 
with initial conditions distributed equally over the stable circle for (JTH2J) . Kicks in the radial direction alone are applied at 
times t = 0, 0.25, 0.5, 0.75. A comparison between a = and a = 2 illustrates the a-induced stretch-and-fold action in the laser 
phase space. 




FIG. 11: Phase-space stretching along the limit cycle shown as time evolution of the phase difference between two trajectories 
starting at different isochrones. The three curves are obtained using (dashed) Eqs. {THIJ with a — 0, (solid) Eqs. (JT}|2jl with 
a = 3, and (dotted) Eq. ((5| with a — 3. Adapted from p3] with permission. 

cycle to allow formation of random strange attractors. Furthermore, we demonstrated that in higher dimensional 
systems, different eigendirections with distinctly different relaxation rates toward the limit cycle could give rise to 
more than one region in the (J, \f2D ext )-]Amie with a random strange attractor. Last but not least, we revealed that 
the enhancement in the instantaneous stretch-and-fold action arising from laser relaxation oscillations results in a 
larger parameter region with a random strange attractor. 

VIII. CONCLUSIONS 

We used the class-B laser model in conjunction with the Landau-Stuart model (Hopf normal form with shear) to 
study noise synchronisation and loss of synchrony via shear-induced stochastic d-bifurcations. 

We defined noise synchronisation in terms of pullback convergence of random attractors and showed that a nonlinear 
oscillator can synchronise to stochastic external forcing. However, the parameter region with synchronous dynamics 
becomes interrupted with a single or multiple intervals of asynchronous dynamics if amplitude-phase coupling or 
shear is sufficiently large. Stability analysis shows that the synchronous solution represented by a random sink 
loses stability via stochastic d-bifurcation to a random strange attractor. We performed a systematic study of this 
bifurcation with dependence on the three parameters: the Hopf bifurcation parameter (laser pump), the amount of 
shear (laser linewidth enhancement factor), and the stochastic forcing strength. 

In this way, we uncovered a vast parameter region with random strange attractors that are induced purely by 
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stochastic forcing. More specifically, in the three-dimensional parameter space, the two-dimensional surface of the 
stochastic bifurcation originates from the half-line of the deterministic Hopf bifurcation. In the plane of stochastic 
forcing strength and Hopf bifurcation parameter, one finds stochastic d-bifurcation curve(s) bounding region(s) of 
random strange attractors if the amount of shear is sufficiently large. The shape of d-bifurcation curves is determined 
by the type and rate of the relaxation toward the limit cycle in the unforced oscillator. Near the Hopf bifurcation 
and provided that stochastic forcing is weak enough, the d-bifurcation curves satisfy the numerically uncovered power 
law (|21[) . However, as the stochastic forcing strength increases, there might be deviations from this law. The deviations 
arise because different oscillators experience different effects of higher-order terms and additional degrees of freedom 
on the relaxation toward the limit cycle. In the laser example, the d-bifurcation curves deviate from the simple power 
law (pH]) so that the region of a random strange attractor splits up and expands toward smaller values of shear as the 
forcing strength increases. We intuitively explained these results by demonstrating that the shear-induced stretch- 
and-fold action in the oscillator's phase space facilitates creation of horseshoes and strange attractors in response 
to external forcing. Furthermore, we showed that the stretch-and-fold action can be greatly enhanced by damped 
relaxation oscillations in the laser model, causing the deviation from the simple power law. 
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